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Abstract 



^ ■ We consider the B-ball cosmology of the MSSM in the context of D-term infia- 

H ; 

, tion models where the reheating temperature is determined by the Affleck-Dine 

mechanism to be of the order of 1 GeV. We show that such a low reheating tem- 
perature can arise quite naturally as the result of a symmetry which is required 
to maintain the flatness of the inflaton potential. In this case the B-balls will 
decay well after the freeze-out temperature of the LSP, allowing baryons and 
cold dark matter to originate primarily from B-ball decays. 
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The MSSM [|I| is known to admit non-topological solitons in its spectrum 0; in 
particular, B-balls (Q-balls P, Q with baryonic charge) can exist. In a cosmological 
scenario which includes inflation these can be copiously produced by the breakdown 
of scalar condensates along flat directions of the MSSM. In the case of gauge-mediated 
SUSY breaking, extensively discussed in references [5-10], the B-balls are stable and 
could account for cold dark matter [Q. In the case of gravity-mediated breaking, 
studied in ITT, 0, the B-balls are unstableQbut nevertheless, if they can survive ther- 



malization, they are typically long-lived enough to decay much after the electroweak 



phase transition, leading to a variant of the Affleck-Dine mechanism [14| which we 



call B-ball Baryogenesis (BBB) |TT|, [T2| . The requirement that they can survive ther- 
malization implies that the B-balls in R-parity conserving models originate from a 
d=6 AD condensate and imposes an upper bound on the reheating temperature of 
10^^^ GeV [|ri|, |12||. Such B-balls can protect a B asymmetry originating in the AD 
condensate from the effects of additional B-L violating interactions, which would oth- 
erwise wash out the B asymmetry when combined with anomalous B-l-L violation [|ll[|. 
In addition, if the reheating temperature is sufficiently low that the B-balls decay be- 
low the freeze-out temperature of the lightest SUSY particle (LSP), then cold dark 
matter can mostly come from B-ball decays rather than from thermal relics. This 
opens up the possibility of relating the number density of dark matter particles to 
that of the baryons, allowing for an explaination of their observed similarity for the 



case of dark matter particles with weak scale masses [|T^ . Since BBB requires low re- 
heating temperatures, it is important to consider whether low reheating temperatures 
are a natural feature of realistic SUSY inflation models. We will show that in D-term 



inflation models ||T^, |I^ very low reheating temperatures are indeed a natural feature. 

In general, SUSY inflation models can be classifled as either F- or D-term inflation 
models, depending on whether the energy density during inflation comes from the F- 



^ Depending on the details of the hidden sector and its interaction with the observable 
sector, it is conceivable that SUSY breaking could be suppressed for large values of the 
condensate field, leading to stable B-balls in the case of gravity-mediated SUSY breaking 
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or D-term contribution to the scalar potential. However, in F-term inflation models, 
H dependent corrections to the mass terms tend to spoil the flatness of the inflaton 
potential ||17| , |l^. As a result, the most favoured candidate for SUSY inflation at 



present is D-term inflation, which has the form of a hybrid inflation scenario [19 



20[]. An important feature of D-term inflation models is that the A-terms receive 



no H dependent corrections [|15|, [T^ (although the mass squared terms do receive 



corrections once inflation ends ||T6[)- This has the important consequence that the 
baryon asymmetry from the d=6 AD condensate will be maximal in a typical domain 
of the Universe, which in turn requires that the reheating temperature is of the order 
of 1 GeV. The main purpose of this paper is to consider whether such low reheating 
temperatures can be a natural feature of realistic D-term inflation models, thus making 
them compatible with BBB and the AD mechanism along d=6 directions. 



The simplest D-term inflation model [0, |T6[ is based on three flelds: a gauge 
singlet 5, which plays the role of the slow-rolling inflaton, and two flelds with opposite 
U{1)fi charges, ip+ and where U{1)fi is the Fayet-Iliopoulos (FI) gauge group. 
The superpotential is given by 

w = xsip+ij- . (1) 

The scalar potential, including the U{1)fi D-term, is then 

V = \X\' + l^^+p + \S^.\') + ^ (l^+P - l^-P + ef , (2) 

where is the FI term and g is the U{1)fi coupling. C, is flxed by cosmic microwave 
background radiation (CMBR) measurements |^ to be 6.6x10^^ GeV [^. The global 



minimum of the potential is at 5* = 0, ip+ = 0, = ^. However, for S > Scru = g(,/X, 
the minimum is at 4'+ = ip- = and there is a non-zero energy density p = 
Although there are no tree-level corrections to the flatness of the S potential due to 
SUSY breaking by the energy density during inflation, there will be an 5* potential 
coming from 1-loop corrections and from higher order terms in in the Kahler 

potential, where M* = Mpi/^/Sn is the supergravity mass scale The 1-loop term, 
due to the splitting in mass of the bosonic components of 'ip± from their fermionic 
partners, is given by AV = a'^C,'^ln{\XS\/ fi) 0, IT^|. 



In this model, slow-rolling of the S scalar fails well before Scrit, at Sf ~ ^JaM^ /2tt. 
In order to account for the observed homogeneity and isotropy of the Universe, the 
slow-roll of 5* must begin at least at 5*55, where 

S,,= (55^ + S?)"^0.9M.. (3) 

For such large values of S, it is likely that higher order superpotential terms will 
cause the 5* scalar potential to deviate from flatness. For example, a term of the form 
AW = nS"^ / M^~^ with S ~ Mt, produces too large a deviation from the observed 
CMBR unless k ^ 10^^. Since we would expect the natural mass scale of the physical 
interactions coming from this superpotential term to be M*, the natural value of k will 
be of the order of 1/m!. Thus in order to have a naturally fiat potential during inflation, 
all superpotential terms with m < 9 must be suppressed by discrete symmetries or R- 
symmetries [|16|. R-symmetries are particularly effective in doing this, since they can 



completely eliminate higher order superpotential terms which are purely a function 
of S as well dangerous gauge kinetic terms of the form S'^W^Wa 0, (There is 
a danger that a global symmetry could be broken by non-perturbative gravitational 
effects ^ , in which case a discrete gauge symmetry would be preferred. However, 
such effects are very sensitive to the details of the gravity theory on small scales, and do 
not generally exclude models with global symmetries [^]). Depending on the charge 
of S and the MSSM fields with respect to the R-symmetry or discrete symmetry, this 
suppression of higher order operators may also suppress the couplings of S to the 
MSSM sector. This, as we will show, can lead to very low reheating temperatures 
after the S oscillations decay. 

In D-term inflation models there are two seperate periods of reheating, correspond- 
ing to the decay of the ip- field (at temperature T^) and of the inflaton S fl^. The 
initial values of the scalar field expectation values at the end of inflation are given 
hy = — ipJ) = ^ and S ~ Scru = fi'^/A. The masses of the scalars ip'_, ip+ 
and S are M^/ = and Ms = M^^ = A^. Thus the initial energy density in the 
oscillating ^p'_ and 5* fields will be essentially the same, ps ~ p^' ~ g^^'^- This has 
the consequence that the Universe will be matter dominated by scalar field oscillations 
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until both condensates have decayed. The immediate ip'- decay leaves the Universe at 
a temperature ~ g^^'^C,- The Universe is subsequently matter dominated, but with 
an initial radiation density much larger than that expected purely from the decays of 
the inflaton S 0. Eventually, at a temperature T^,, the radiation coming from the decay 
of the S inflatons, which has an energy density given by 

'iTsPs 



5H 



(4) 



where Ts is the S decay rate, comes to dominate that remaining from the first stage 
of reheating, so that the subsequent evolution is the same as for the minimal single 
inflaton model of Refs.[|lT], O. The expansion rate at T^,, if*, is given by 



B-balls typically begin to form at Hi ^ Sli^l GeV, where \K\ is the coefficient of 
the logarithmic radiative correction to the B-ball scalar potential; typically \K\ ^ 



0.01 — 0.1 pTI , |T2|. Thus H^, is typically smaller than Hi, especially for the case of 
small T^. This means that the temperature of the Universe when the AD condensate 



forms (at H ^ 100 GeV 0, O) and when the B-balls begin to form at Hi will be 



somewhat higher than in the single inflaton model. However, there is no danger of 
the B-balls being thermalized so long as the AD condensate from which they form is 
not itself thermalized. This is because the background of non-relativistic squarks and 
anti-squarks which exists after the AD condensate collapses to B-balls will give large 
effective masses T) to the particles which could otherwise thermalize the B-balls 
[P^. Therefore the formation and evolution of the B-balls in the D-term inflation 
scenario will be effectively the same as in the minimal single inflaton model. 

We next consider the decay of the S scalar oscillations. As noted above, the self- 
couplings of S in the superpotential must be highly suppressed by a symmetry in order 



^The S condensate could be thermalized by this radiation density, via scattering with 
light MSSM scalars once the U{1)fi gauge fields are integrated out |l^. However, thanks 
to the strong dependence of the scattering rate on A {Tscatt oc A^), this will not occur so long 
as A < 0.5. 
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to ensure that the scalar potential is flat enough for successful inflation. This in turn 
can suppress the couplings of S to MSSM fields. Suppose that the lowest dimension 
superpotential coupling of the S scalars to the MSSM fields is of the form 



W = ^ (6) 



where represents the MSSM fields. This gives, for example, a decay rate to r 
particles of the form 



r, « (^^j .^AM. , (7) 
where Pr represents the phase space factor. The decay temperature is then 



where during radiation domination H = kxT'^/Mpi with ~ 17. With Ms = XC, this 
gives 

T| ^ 4 X 10^5(2.8 X 10-''y~''/W~^/^(3y^ . (9) 

Low reheating temperatures can arise as a result of a "mismatch" between the trans- 
formation properties of the fields of the infiaton sector {S, ip^ and ip^) and the MSSM 
sector under the symmetry. As a particularly simple example, consider an R-symmetry 
under which the infiaton sector fields have charges which are multiples of 1/2 and the 
MSSM sector have charges which are multiples of 1/3. An example which allows the 
XSip+ip- superpotential term (allowed terms have a total R-charge equal to 2 |jl|) 
but eliminates the higher order S superpotential terms is defined by R{S) = — n, 
R{iIj+) = R{ip-) = (n + 2)/2 for the infiaton sector fields and R{Q) = R{L) = 5/3, 
R{u') = R{d^) = R{e^) = 1/3 and R{Hu) = R{Hd) = for the MSSM fields. (The 
MSSM charges have been chosen in order to allow the higher order terms which stabi- 
lize the d=6 u'^d'^d'^ D-flat direction). In Table 1 we give the lowest order R- invariant 
operators coupling S to the MSSM fields for the case where n is an positive integer 
with n < 3. (For > 3 the operators are of higher order). 
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Table 1. R-invariant operators. 



n 


r 




1 


6 




2 


4 


QQu'd' 


3 


6 





For example, with n = 1 the lowest order superpotential coupling of 5* to the MSSM 
fields has r = 6. This gives a reheating temperature 

9/2 



150 [^j Kf3l/'GeV 



(10) 



where /?6 ~ 10^ . Thus a low reheating temperature, as low as 1 GeV or less, is a 
natural feature of this model. 

In general, we should also consider corrections to the Kahler potential and the 
gauge kinetic term which define the full supergravity theory It is straightforward 
to show that these non-minimal corrections do not introduce any operators linear in S 
of lower order in M~^. It is also possible that such non- minimal corrections could allow 



for rapid decay via parametric resonance |27, 28|. However, it is straightforward to 
show 1^ that the resulting resonant decay of the infiaton will be completely negligible 
so long as X/g^^"^ ^ 0.03. For values of A and g not too small compared with 1 this 
condition will be easily satisfied. 

We next consider the need for low reheating temperatures in d=6 AD baryogenesis. 
The MSSM has many D-fiat directions in the scalar potential, along which the fiatness 
is lifted only by SUSY breaking terms and non-renormalizible superpotential terms 
[E^]. The potential of the complex AD scalar $ will have the form 



f/($) 



cH^) 1 + is: log 



|<|)|2 



where d is the dimension of the non-renormalizible term in the superpotential and 
cH^ gives the H correction to the scalar mass (with c positive and of the order of 
1). The logarithmic correction to the scalar mass term, which is expected in SUSY 
models, is crucial for the growth of perturbations of the AD field and the formation 



of B-balls, which occurs if < 0. Of particular interest to us here will be the 
u'^d'^d'^ direction, for which < and along which a d=6 superpotential term of the 
form {u'^d'^d'^Y is allowed by the R-symmetry of the previously discussed model. This 
superpotential term (via the associated A-term in the scalar potential) is responsible 
for introducing the B- and CP-violation necessary to produce the baryon asymmetry 
via the AD mechanism At cH^ > m^, the scalar field will be at the minimum of its 
potential, with an expectation value given by |$| = $o ^ {cH'^M^ /5\^Y^^. However, 
in the absence of order H corrections to the A-term, the phase of $ will be completely 
undetermined. Once cH^ ^ m^, the AD field will start to oscillate. The A-term then 
creates a phase shift 6 of the order of 1 between the real and imaginary parts of the 
oscillating AD field. With initially $ = $06*'''^^, where 6cp is the random initial 
phase determined during inflation, the initial B asymmetry in the AD condensate at 
Ho ~ mjc^l'^ is given by n^o = em^^/S, where e = Sin25cpSin5. This gives for the 
baryon to entropy ratio of the d=6 AD condensate, 

ne 27mBoTi /0.008N '/'^ /lOO GeV\i'2 / r| \ 

'te-T = ^/51vJr m l^^J [WGevj' 

where we have assumed that A ~ 1/5!, such that the strength of the physical self- 
interactions of $ is determined by the mass scale M^,, and we have taken 6 ~ 6cp ~ 1. 
The corresponding reheating temperature for the d=6 AD condensate is 

Nucleosynthesis implies that ^ (3 — 8) x 10^^^ for an Q = 1 Universe. Thus, 
for 6cp ~ 1, we find that ^ 1/c GeV, which can reasonably take values in the 
range O.lGeV to lOGeV, depending on c, with an "average" value of around 1 GeV. 
This is consistent with the range of values attainable in R-symmetric D-term inflation 
models. 

A fraction of the total B asymmetry, fs, expected to be not very small compared 



with 1 , ends up in the form B-balls once the AD condensate collapses to a mixture 
of B-balls and non-relativistic squarks ai H ^ Hi. The B-balls have very large charges, 
B ^ 10^'^/b(1 GeV/T^), and, for ^ 1 GeV, decay at a temperature in the range 



1 MeV to 1 GeV |T2|. In this case the B-balls will decay below the LSP freeze-out 
temperature, Tf ^ m^/20, where m^^ is the LSP neutralino mass and the LSP 
density is likely to be dominated by LSPs coming from B-ball decays [1^. This is 
particularly true in the present model in which there is a low reheating temperature, 
since, if < Tj, the thermal relic LSP density will be strongly suppressed (by 
an additional factor {T^/TfY) by the period of matter domination between Tf and 
T^. The domination of the LSP density by B-ball decays is, in itself, an important 
consequence of B-balls, requiring a new analysis of dark matter constraints on the 
MSSM ||3^. It also makes it possible to account for the similarity of the number 
densities of baryons and LSPs, since, if fs is not very small compared with 1, the 
number density of LSPs from B-ball decays will be naturally similar to that of baryons 

Thus, assuming nothing more than a low energy sector consisting of the MSSM 
fields together with a period of primordial D-term inflation, we can, in principle, nat- 
urally account for all the main observation of cosmology; the baryon asymmetry (via 
the Affleck-Dine condensate), cold dark matter and the baryon to dark matter number 
ratio (LSP neutralinos from B-ball decays) and the homogeneity of the Universe to- 
gether with the CMBR temperature fluctuations (D-term inflation). This provides us 
with an economical and effective model for SUSY cosmology. The resulting scenario 
radically departs from the conventional radiation dominated view of SUSY cosmology 
at low temperatures, with the electroweak phase transition playing no signiflcant role 
in baryogenesis and late decaying B-balls opening up new possibilities for SUSY dark 
matter. 

Acknowledgements 

This work has been supported by the Academy of Finland and by a European Union 
Marie Curie Fellowship under EU contract number ERBFM-BICT960567. 



8 



References 

[1] H.P.Nilles, Phys. Rep. 110 (1984) 1. 
[2] A.Kusenko, Phys. Lett. 405 (1997) 108. 
[3] S.Coleman, Nucl. Phys. B262 (1985) 263. 

[4] A.Cohen, S.Coleman, H.Georgi and A.Manohar, Nucl. Phys. B272 (1986) 301. 

[5] A.Kusenko, Phys. Lett. B404 (1997) 285. 

[6] A.Kusenko and M.Shaposhnikov, Phys. Lett. B418 (1998) 46. 

[7] A.Kusenko, V.Kuzmin, M.Shaposhnikov and P.G.Tinyakov, Phys. Rev. Lett. 80 
(1998) 3185. 

A.Kusenko, M.Shaposhnikov, P.G.Tinyakov and I.G.Tkachev, Phys. Lett. B423 
(1998) 104. 

[9] G.Dvali, A.Kusenko and M.Shaposhnikov, Phys. Lett. B417 (1998) 99. 
[10] M.Laine and M.Shaposhnikov, \hep-ph/ 980423^ . 
[11] K.Enqvist and J.McDonald, Phys. Lett. B425 (1998) 309. 
[12] K.Enqvist and J.McDonald, \hep-ph/980338dl . 



[13] A.Kusenko, private communication. 

[14] I.A.Affleck and M.Dine, Nucl. Phys. B249 (1985) 361. 

[15] E.Halyo, Phys. Lett. B387 (1996) 43; P.Binetruy and G.Dvah, Phys. Lett. B388 
(1996) 241. 



[16] C.Kolda and J.March-Russell, \hep-ph/9802358 



[17] M.Dine, W.Fischler and D.Nemeschansky, Phys. Lett. B136 (1984) 169; 
G.D.Coughlan, R.Holman, P.Ramond and G.G.Ross, Phys. Lett. B140 (1984) 
44. 



[18] E.Copeland, A.Liddle, D.Lyth, E.Stewart and D.Wands, Phys. Rev. D49 (1994) 
6410; E.D.Stewart, Phys. Rev. D51 (1995) 6847; M.Dine, L.Randall and 
S.Thomas, Phys. Rev. Lett. 75 (1995) 398; G.Dvali, fiep-ph/950325^ . 



[19] K.Enqvist and D.V.Nanopoulos, Phys. Lett. B142 (1984) 349. 
[20] A.Linde, Phys. Rev. D49 (1994) 748. 

[21] G.F.Smoot et al., Ap. J. 396 (1992) LI; A.Stebbins, {astro-ph/gjOWTS . 
[22] D.Lyth and A.Riotto, Phys. Lett. B412 (1997) 28. 
[23] D.Lyth, Phys. Lett. B419 (1998) 57. 

[24] S.Giddings and A.Strominger, Nud. Phys. B307 (1988) 854; S.Coleman, Nucl. 
Phys. B310 (1988) 643; G.Gilbert Nucl. Phys. B328 (1988) 159; S.Coleman and 
K.Lee, Phys. Lett. B221 (1989) 242; Nucl. Phys. B329 (1990) 387; L.Krauss and 
F.Wilczek, Phys. Rev. Lett. 62 (1989) 1221. 

[25] M.Kamionkowski and J.March-Russell, Phys. Lett. B282 (1992) 137; R.Holman, 
S. D.H.Hsu, T.W.Kephart, E.W.Kolb, R.Watkins and L.M.Widrow, Phys. Lett. 
B282 (1992) 132. 

[26] R.Kallosh, A.Linde, D.Linde and L.Susskind, Phys. Rev. D52 (1995) 912. 

[27] L.Kofman, A.D.Linde and A.A.Starobinsky, Phys. Rev. Lett. 73 (1994) 3195; 
D.T.Son, Phys. Rev. D54 (1996) 3745. 

[28] T.Prokopec and T.G.Roos, Phys. Rev. D55 (1997) 3768; I.Zlatev, G.Huey and 
RJ.Steinhardt, Phys. Rev. D57 (1998) 2152. 



M.Dine, L.Randall and S.Thomas, Nucl. Phys. B458 (1996) 291. 
[30] S.Sarkar, Rep. Prog. Phys. 59 (1996) 1493. 

[31] J.Ellis, J.S.Hagelin, D.V.Nanopoulos, K.Olive and M.Srednicki, Nucl. Phys. B238 
(1984) 453; G.Jungman, M.Kamionkowski and K.Griest, Phys. Rep. 267 (1996) 
195. 

10 



K.Enqvist and J.McDonald, ]hep-ph/98U7M^ . 



11 



